ICRR-Report-317-94-12 
April 1994 



Flavor mixing in the gluino coupling 
and the nucleon decay 

Tom Goto, Takeshi Nihei and Jiro Arafune 



Institute for Cosmic Ray Research, University of Tokyo, 
Midori-cho, Tanashi-shi, Tokyo 188 JAPAN 



Abstract 

Flavor mixing in the quark-squark-gluino coupling is studied for the 
minimal SU(5) SUGRA-GUT model and applied to evaluation of the nu- 
cleon lifetime. All off-diagonal (generation mixing) elements of Yukawa 
coupling matrices and of squark/slepton mass matrices are included in 
solving numerically one-loop renormalization group equations for MSSM 
parameters, and the parameter region consistent with the radiative elec- 
troweak symmetry breaking condition is searched. It is shown that the 
flavor mixing in the gluino coupling for a large tan (3 is of the same order 
of magnitude as the corresponding Kobayashi-Maskawa matrix element 
in both up-type and down-type sector. There exist parameter regions 
where the nucleon decay amplitudes for charged lepton modes are dom- 
inated by the gluino dressing process, while for all the examined regions 
the neutrino mode amplitudes are dominated by the wino dressing over 
the gluino dressing. 



I Introduction 



SU(5) supersymmetric (SUSY) grand unified theory (GUT) is an attractive can- 
didate for the unified theory of strong and electroweak interactions. The analyses 
of the gauge coupling unification |J suggest the validity of the minimal supersym- 
metric standard model (MSSM) just above the electroweak scale ~ 100 GeV and 
the unification of SU(3)xSU(2)xU(l) gauge group into a simple SU(5) at the GUT 
scale M x ~ 10 16 GeV . 

One of the features of MSSM/SUSY-GUT is the existence of soft SUSY break- 
ing. It gives quarks (or leptons) and their superpartners different mass matrices in 
the generation (flavor) space. In results, due to the discrepancy between the mass 
diagonalizing bases of quarks and those of squarks, a generation mixing occurs in 
the quark-squark-gaugino coupling (gluino coupling, in particular), which may give 
considerable contributions to the nucleon decay, flavor changing neutral currents 
(FCNC) and other various phenomena in SUSY-GUT. The flavor mixing in the 
gaugino coupling plays an important role in the nucleon decay, since its amplitude 
is dominated by the dimension five interaction followed by the gaugino "dressing" 
process fU in the minimal SU(5) SUSY-GUT model. However, only a simplified 
treatment is made in the previous analyses of the nucleon decay || [3], |, where 
the diagonal (in the generation space) gluino coupling is assumed leading to the 
negligible contribution from the gluino dressing process. 

The purpose of the present paper is to study the flavor mixing in the gaugino 
coupling extensively and evaluate the contribution of the gluino dressing process to 
the nucleon decay amplitude. Based on the minimal SU(5) supergravity (SUGRA-) 
GUT |7j , we assume that the soft SUSY breaking parameters are "universal" at the 
GUT scale [|[. We include all off-diagonal elements of Yukawa coupling matrices 
and of squark/slepton mass matrices in solving numerically one-loop renormaliza- 
tion group equations (RGEs)Q for all MSSM parameters |10|, |11[] with the universal 
boundary conditions. We then evaluate the effective potential for the Higgs fields at 
the electroweak scale to find a consistent SU(2)xU(l) breaking minimum in accor- 
dance with the radiative electroweak symmetry breaking scenario PD[. The obtained 



mass matrices of all particles are diagonalized to evaluate the flavor mixing in the 

* Two- loop RGEs for soft SUSY breaking parameters are obtained recently which will be 
important for more accurate analysis. 
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gaugino couplings. The mass spectrum and the mixing are then used to calculate the 
nucleon decay amplitudes for various decay modes. It is found that the flavor mixing 
in the gluino coupling depends on tan (3 (ratio of vacuum expectation values of the 
Higgs doublets) and is roughly of the same order of magnitude as the corresponding 
Kobayashi-Maskawa matrix element. As for the nucleon decay, we find the gluino 
dressing process dominates the amplitude for the decay modes containing a charged 
lepton (and a meson) if tan j3 is large and if the gluino mass is much smaller than 
the squark masses. Note that the flavor mixing in the gaugino couplings is studied 
previously in a systematic analysis of FCNC ]ll with a semi-analytic solution of 
the RGEs for small tan/3 (top Yukawa coupling 3> bottom Yukawa coupling). On 
the contrary, our numerical method cover the whole range of tan /3, since all Yukawa 
couplings are taken into account. 

The remaining part of this paper is organized as follows. After introducing 
in the next section the minimal SU(5) SUGRA-GUT model which we consider, 
we formulate the flavor mixing in the gaugino couplings in Sec. |ITT[ In Sec. [IV] the 
nucleon decay amplitudes are obtained with a careful treatment of the flavor mixing. 
The outline of the numerical calculation and the results are presented in Sec. |VJ and 
our conclusions are summarized in Sec. |Vl. 



II Minimal SU(5) SUGRA-GUT model 

Minimal SU(5) SUSY-GUT model contains three generations of matter multiplets 
with 10 and 5 representations, tyf B and respectively, where suffices A, B = 
1,2, ••-,5 are SU(5) indices and i,j = 1,2,3 are the generation labels, and three 
kinds of Higgs multiplets with 5, 5 and 24 representations, H$, H^ A and re- 
spectively. SU(5) and R-parity invariant superpotential Wgut at the GUT scale is 
written as 

W GUT (M X ) = r^t B <b 3A H ZB + ^g l ^ABCDE^t B ^ D Hi 

+XH- 5A (E£ + M8*) Hi + W s (£) . (2.1) 

Here, ^abcde is the totally antisymmetric constant, f l \ = g^ 1 and A are dimen- 
sionless couplings, M is a mass parameter and Wy, is a self-interaction superpotential 
for E^. SU(5) symmetry is spontaneously broken down to SU(3) x SU(2) x U(l) 
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with the nonvanishing vacuum expectation value of the adjoint Higgs Below 
the GUT scale, the model is reduced to MSSM with effective higher dimensional op- 
erators, which are obtained by integrating out the superheavy particles in Eq. (|2.f|). 
The effective superpotential is then written as 

W cS = Wigsm + W 5 + 0(Mx 2 ) , 

W5 = —— \ -C % [ kl t a bctafiQ°k $1 QT 'Lj-y + C]l kl e abc Ui a Dj C EkUib 
Mo I 2. 

+ (baryon number/lepton number conserving terms) | . (2.2) 

Here, Q, U and E are chiral superfields which contain left-handed quark doublet, 
right-handed up-type quark and right-handed charged lepton respectively, and are 
embedded in \1/ (to be specified in Eq. (|3.4Q ); D and L, which are embedded in <3>, 
contain right-handed down-type quark and left-handed lepton doublet respectively; 
Hi and H 2 are Higgs doublets embedded in if 5 and if 5, respectively. The suffices 
a,b,c = 1, 2, 3 are SU(3) indices and a, [3 = 1,2 are SU(2) indices. M c is the colored 
Higgs mass which is assumed to be O(Mx), while the supersymmetric mass of Higgs 
doublet fi is of the order of the Z boson mass mz- This discrepancy is owing to a 
tree level fine-tuning in the GUT superpotential. At the GUT scale, fn and fi are 
unified and Cl and Cr are written in terms of the Yukawa coupling constants (see 
Sec. |TV|) . Baryon number (and lepton number) violating terms in W5 give dominant 
contributions to the nucleon decay in this model. 

In addition to the supersymmetric Lagrangian to be derived from ( |2.2j ), the 
following soft SUSY breaking terms are included: 

- £ so ft = ( m Q^ tf l Qj + ( m r>)^ d^dj + (m^, u h Uj 

+ (miy, 1% + (m 2 E y. e% + A^i + ^\h\h 2 

+ ^A^qidjhi + AYeJjhi + A^qiU^ - B/j,hih 2 + h. c.} 

+-{M 1 BB + M 2 WW + M 3 GG + h. c.} , (2.3) 

where cji, di, Ui, e«, k, hi and h 2 are scalar components of Qi, Di, Ui, E iy Li, Hi and 
H 2 , respectively, and B, W and G are U(l), SU(2) and SU(3) gauge fermion fields 
(bino, wino and gluino), respectively. SU(2) and SU(3) suffices are omitted in ( |2.3|) 
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for simplicity. We assume that the soft SUSY breaking parameters satisfy simple 
relations at the GUT scale: 



A\ = A 2 . = ml , 

a d = fr)x A xm , A l [ = f L 3 x A x m , A J J = g$ x A x m Q , 

■U; M 2 M 3 - M gX , (2.4) 

where the suffix "X" stands for the value at the GUT scale. The boundary conditions 
( p.4| ) are due to the minimal SUGRA model, where local SUSY is spontaneously 
broken in the hidden sector which couples to the observable sector (SUSY-GUT, 
in the present case) only gravitationally, and hence universal soft SUSY breaking 
terms are induced in the observable sector ||. 

Below the GUT scale, radiative corrections modify all parameters in the su- 
perpotential ( |2.2| ) and the soft SUSY breaking terms (^.3|), as well as three gauge 
coupling constants g\, gi and g% for U(l), SU(2) and SU(3), respectively. The evo- 
lution of the parameters are described by the RGEs [|TT|] . According to the radiative 
SU(2) x U(l) breaking scenario [10] , we numerically solve the RGEs down to the 



electroweak scale rriz and evaluate the effective potential for the neutral Higgs fields: 

U(Higgs) = V + V x , 

V = (fi 2 + A 2 1 )\h 1 \ 2 + (fi 2 + Al)\h 2 \ 2 -(Bfih 1 h 2 + h.c.) 

92l+gh \h 1 \ 2 -\h 2 \ 2 ^ 



where Str means the supertrace and M. includes all (s)quark and (s)lepton masses. 
Then the electroweak symmetry breaking condition 



(hi) = vcos/3, (h 2 ) = v sin/3 , (2.6) 

= 92 v 2 , 
2 cos 2 6 W ' 



rrir 



is imposed. 
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Ill Flavor mixing in the gluino coupling 



In order to discuss the flavor mixing in the gluino coupling, we have to diagonalize 
the mass matrices for quarks and squarks. Throughout the calculation hereafter, 
we choose the basis in the generation space for the superfields so that the Yukawa 
coupling constants for up-type quarks and leptons should be diagonalized at the 
electroweak scale. The Yukawa terms in (|2.2| ) are then written as 

^Yukawa(mz) = f% (^m)J QiDjH x + fj^E^H, + $QiU 3 H 2 , (3.1) 

where the notation " " " stands for a diagonal matrix and Vkm is the Kobayashi- 
Maskawa matrix. All eigenvalues of /d, fi and gu are taken to be real positive. 
Since this choice of the basis is different from that in the GUT superpotential ( [2.1D , 
a re-diagonalization of the Yukawa couplings at the GUT scale is needed in order 
to find the unification condition of Jd and fi and the relation between the Yukawa 
coupling constants and the dimension-five coupling constants Cl,r- JUyukawa at the 
GUT scale is diagonalized with appropriate unitary matrices Uq \ Ud, Ue, Ul, Uq^ 
and Uxj\ 

^Yukawa(M x ) = f^xQiDjHx + f^ x E i L j H 1 + gij x QiUjH 2 (3.2a) 
= fox {U^)l (U D )/ QtPjHx + f k L l x (U B )* (U^ E t L d H, 
+g k u l x (u i ^)*(U u )/Q l U J H 2 . (3.2b) 

The unification condition of fjj and ft, is then written as[] 

fix = fix • (3-3) 
The matter multiplets are accommodated into the SU(5) multiptets as 

<= {Qi, (u^ptuu)'^, (yWuB^Ej} , 

*i <= {A, (U D U L )'L^ , (3.4) 

^The "unification" of fox and /lx, with quark masses given in Ref. [|l2|, however, is not so 
satisfactory numerically in the first and the second generations as the gauge coupling unification. 
We ignore the difference in the present calculation, since we may still have ambiguities of the 
renormalization effect in the very low energy region. 
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and the GUT Yukawa coupling constants in Eq. ( |2.1| ) are expressed with those in 
Eq. ( |3T2il ) as 



J — Jr 



DX , 



g' J 



9b k x(u}jPU^)^ , (3.5) 

where P is a diagonal phase matrix which cannot be absorbed by field redefinitions 
in the colored Higgs coupling [|13| . 

The origin of the flavor mixing in the gluino coupling lies in the difference 
between the mass basis for quarks and that for squarks. The mass matrix for up- 
type squarks is expressed as 

£(s-up mass) = {q u , v))M\ [ % 




= (MuMp) 1 . + (m 2 Q y. + m 2 w cos 2/3 (- - - tan 2 6w) 6 



(m 2 ^ = (MtM c/ )' + (m 2 / )'+m 2 w cos2/?(^-tan 2 ^j^ , 

( m iij) 3 = ^Mjj cot (3 + A^v sin (3 , 

mRL = m\ R , (3.6) 

where M\j is the up-type quark mass matrix M{j = g}}v sin (3 and q u is the up-type 
component of the SU(2) doublet q. The squark mass matrix A4 U is not diagonalized 
with the quark mass basis (|3.1|) since off-diagonal elements are induced in the soft 
SUSY breaking parameter matrices due to the renormalization effect. Squark mass 
basis is obtained by diagonalizing (|3.6|) with a 6x6 unitary matrix Ujj: 



u> 



rt \si2i 



(UuYjUj, I = 1, 2, 6 




UjjMiUu = diagonal , (3.7) 



where u\ is the mass eigenstate of up-type squark. We define the numbering of u\ 
such that the mixing of uj is the largest in vlj. Accordingly we call u[, u' 2 , ■ ■ ■, u' 6 
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as ul-, Cl, £l, ur, cr and Ir respectively in the later discussions. The mass bases 
of down- type squarks and charged sleptons are obtained in the same way with 6x6 
unitary matrices Ud arid Ue, respectively. Notice that no generation mixing occurs 
in the lepton/slepton sector since the right-handed (s)neutrino does not exist in the 
minimal model; the nonvanishing off-diagonal elements of the slepton mass matrix 
are left-right mixing components only. Consequently, quark-squark-gluino coupling 
is written as 

£ int (gluino) = -iV2g 3 (d^ I ((u D ) k i (V KM )iGd Lj +^ 

+ u'V (^Uv) j I Gu Lj + (UuY^GurM +h. c. 
= -iV2g 3 {d'V (jfj 1 ^ Gd Lj + (Uj,)^ Gd Rj >j 

+ «' f/ (jUuYjGuLj + (UuY^GuR^+h. c. , (3.8) 

with the definition of U' D as 

uii and du in ( p.8[) are left-handed quarks of mass eigenstates which compose the 
SU(2) doublet as 

and uri and dRi are the fermion components of C/j and Di, respectively. Similar 
flavor mixing formulae are obtained for other gaugino (wino and bino) coupling 
terms. 

IV Nucleon decay with dimension five operators 

As mentioned in Sec. [HJ the nucleon decay amplitude in the minimal SU(5) SUSY- 
GUT model is dominated by the dimension five operators P| induced by colored 
higgsino/Higgs exchanges. Since the dimension five operators are made from two 
fermion (quark/lepton) and two boson (squark/slepton) component fields, effective 
baryon number violating four-fermion operators are generated by one loop "dress- 
ing" diagrams which involve gauginos or higgsinos (see Fig. |l|). In the present 
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calculation, only the gluino dressing and the charged wino dressing diagrams are 
included; contributions from higgsino dressing diagrams are negligibly small due to 
the small Yukawa couplings of light quarks (u, d, s), compared to the SU(2) gauge 
coupling g2] neutral wino and bino coupling have the same flavor mixing structure 
as that in the gluino coupling, hence their contributions are smaller than that from 
the gluino dressing. 

The dimension five coupling constants Cl and Cr of (|2.2|) at the GUT scale 
are written in terms of the Yukawa coupling constants (see (|3.4j ) and (|3.5|) ): 

Cg? = fox {ulU L y m g$ x (UIPU^)[ , 

ess = /^K^^L^K"^)! ■ (4- 1 ) 

Note that this relation with the index "X" removed does not hold true at the 
electroweak scale. The effective baryon number violating four-fermion operators at 
the electroweak scale are written as 

C ea (AB = ±1) = (ci jk \G) + C^ k \W))e a U<Ai)(d c Li u Lj ) 

+ (Cf\G) + &J k \W)) e a6c K*4)K^) 

+ (right-handed quark/lepton) + h. c. , (4.2) 



where C UjB are calculated as follows with use of the numerical values of Cl and Cr 
at the electroweak scale to be obtained through their RGEs, squark/slepton mass 
eigenvalues and the mixing matrices in the gaugino couplings"'': 




* Contributions from Crs and higgsino/neutralino dressings are estimated to be small and 
neglected in the present calculations. 
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crm = -§-{(ct mn - cr m ) (yhY m {u ] D ) J n F w {u'j,d'j) (u Djj 



+ (cr H - Cf m ) (U^F^, v\) (U D ))} . (4.3) 

Here, (U' D )) is defined in ( |3T9|) and (U^)} = (Uu)\(Vkm))- Fq and are obtained 
by the loop integral || f|. |5[ : 

Fc(kh) = F(m h ,m f ~ 2 ;M 3 ) , 

F W {kh) = {U^Fim^m^MDiui) 1 ^ (4.4) 

~ , , , s 1 M ( m? . m? . m| \ 

I67H mf — \mf — M J m^ — M 2 M 2 J 

where Z7_, C/+ are 2x2 unitary matrices which diagonalize the chargino mass matrix 
and M± (a = 1,2) are its eigenvalues: 

/ M 2 V2m w smf3 \ 

M chargino = I / - 

V — V2mvy cos p — /i y 



The low energy QCD correction between and 1 GeV is taken into account 
in order to evaluate the four fermion operators in the next section. The quark 
Lagrangian at ~ 1 GeV is then converted to the hadron chiral Lagrangian with 
AB = ±1 terms P, [PI with use of the matrix element 



0\e abc (d a L u b L )u c L \p) = (3 P N L } (4.6) 

where Nl is a left-handed proton wave function; it enables us to evaluate the partial 
lifetimes of the nucleon decay. 

V Numerical calculations 

According to the framework described in the previous sections, we calculate the 
flavor (and left-right) mixing in gaugino couplings and the nucleon partial lifetimes 
in a five- dimensional parameter space {m t0 p, tan/5, m , M g x, Ax}, where a di- 
mensioned A parameter is defined as Ax = Axttiq. Actual calculations are made 
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in the following procedure. At first, m top and tan/3 (at the electroweak scale) are 
fixed. Using the numerical values of light quark masses and the Kobayashi-M askawa 



mixing angles given in literatures [T2|, [L5| with the above fixed m top and tan/3, we 
evaluate the Yukawa coupling constants at the electroweak scale Q3.1|) . QCD cor- 
rections below the electroweak scale for quark masses other than m top are included 
at the one-loop level. Next, the RGEs for the dimensionless parameters i.e. , the 
gauge coupling constants and the Yukawa coupling constants are solved upward to 
the GUT scale with the boundary conditions at the electroweak scale. At the GUT 
scale, the Yukawa coupling constants are re-diagonalized to obtain the boundary con- 
ditions for the dimension- five coupling constants (|4.1| ). Then the RGEs for the soft 
SUSY breaking parameters and dimension-five coupling constants are solved down- 



ward with the boundary conditions ( |2.4| ) and (|4.1|) . Since the RGEs are linear for 



the dimensionful parameters, all soft SUSY breaking parameters at the electroweak 
scale are written as linear combinations of the initial parameters (m , M gX , A x ) 

E- 

rh](m z ) = cuml + c 2 iM gX + c 3I A 2 x + c AI M gX A x , 

Mj(m z ) = dxjM gX + d 2J A x , (5.1) 

where rhj and Mj are collective notations for the soft SUSY breaking parameters 
of mass dimension two (ttiodulej ^12) an d one (-^1,2,3, A DtUjL ), respectively. 
The coefficients c's and cfs are implicit functions of the gauge couplings and the 
Yukawa couplings and are determined numerically by solving the RGEs with four 
cases of boundary conditions (m , M gX , A x ) = (1, 0, 0), (0, 1, 0), (0, 0, 1) 
and (0, 1, 1). Once the coefficients are obtained, the values of soft SUSY breaking 
parameters at the electroweak scale for given (m , M gX , A x ) are evaluated with 
the formulae ( |5.1D , and it is easy to scan the three-dimensional parameter space 
{m , M gX , A x } for fixed m top and tan/? with this method. The next step is to 
evaluate the remaining two parameters /i and B with the electroelectroweak SU(2) 
x U(l) symmetry breaking condition. The requirement that the minimum of the 
Higgs potential (|2.5|) gives the vacuum expectation values (|2.(j|) leads to 



A§-A? _ A? + A| _ 1 

^ 2 cos 2(3 2 2 z 
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cos p — — — - sm p 



(5.2) 

h\=v cos p 
h,2=v sin P 



v cos 2/3 \dh\ dh 2 
Bii = ^ + A? + A|\ sin2/3 

+- f^isin^ + |^cos^ . (5.3) 

f a/12 / /ii=-ucos,3 

li2=v sin /3 

Notice that the solution of the equation ( |5.2| ) for /i cannot be written in a simple 
formula since the one-loop part of the Higgs potential V% depends on /x. We solve 
( |5.2[ ) numerically for both signs of /i and then calculate B with ( |5.3| ). Since all 
the MSSM parameters and the dimension-five coupling constants at the electroweak 
scale for a given parameter set (m top , tan/3, m , M 9 x, Ax-) ar e thus determined, 
the mass spectrum of all superparticles and the mixing matrices in the gaugino 
couplings are obtained by diagonalizing their mass matrices. 

We investigate the parameter space {mo, M g x, Ax} within the range 10 GeV < 
m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5 for each combina- 
tion of m top = 120, 150 or 180 GeV and tan/3 = 2, 10, 30 or 50. For m top = 180 
GeV and tan/3 = 2, the top Yukawa coupling diverges below the GUT scale when 
solving the RGEs. For tan (3 = 50, no consistent radiative breaking solution is found 
for any m t0 p. Figs. ^a| - 12| are histograms for the specified off-diagonal elements of 
the gluino coupling matrices Uu and U' D for m top = 150 GeV. The magnitudes of 
the generation mixing in the right-right and left-right sectors are small compared to 
the corresponding left-left elements. For a small tan /3, the left-left elements of U' D 
are approximately equal to the Kobayashi-Maskawa matrix elements: {U' D )\ ~ V c d, 
(U' D )l ~ V t d, etc. in most of the parameter space, while the corresponding off- 
diagonal elements of Uu are small: the mass matrices of up-type quarks, up-type 
squarks and down-type squarks are diagonalized in the same basis. This agrees with 



the conclusion of Ref. |IT| where those mixing matrices are semi-analytically ob- 
tained with an assumption that the bottom Yukawa coupling is much smaller than 
the top Yukawa coupling, which is applicable for small tan/3. On the other hand, 
for a large tan/3, we find that nonvanishing off-diagonal elements in the left-left 
sector of Uu arises with the same order of magnitudes as Vkm, which contributes 
significantly to the charged lepton decay modes (see below). Off-diagonal elements 
of U' D are smaller than those for small tan (3. The qualitative behavior of the mixing 
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matrices is rather independent of the different values of m top - 

Using the obtained values of superpartner masses and mixing matrices, we 
evaluate C's in ( |4.2|) with the formula ( |4.3| ). We then take into account of the low 
energy QCD correction to C's at the one-loop level|. We take the chiral Lagrangian 
factors given in Ref. to derive amplitudes of various decay modes from (|4.2j) . A 
large uncertainty of the nucleon lifetime comes from the numerical values of (3 P and 



Mc- ftp is calculated with various methods [0, which give 

0.003GeV 3 <P P < 0.03GeV 3 , 

and Ref. §, |8| shows 

2.2 x 10 13 GeV < M c < 2.3 x 10 17 GeV . 

Here, we take a small value of f3 p = 0.003 GeV 3 and a large value of Mc = 10 17 
GeV so that we have a longer nucleon lifetime for the safety of the later arguments. 
Fig. [5] shows the partial lifetime for each nucleon decay mode with fixed m top = 150 
GeV and tan (3 = 2. The range of each lifetime comes mainly from the ranges of the 
soft SUSY breaking parameters. As can be seen in the figure, Kv decay modes are 
dominant and most severely constrained by the experiments [|19|, [2(J as 
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Our main concern is to study the contribution of the gluino dressing diagrams. 
To do that, we compare r(wino) with r(total), where r(wino) is a partial lifetime 
calculated with only the wino dressing diagrams taken into account and r(total) 
is that calculated with both wino and gluino dressing diagrams combined. The 
results for m top = 150 GeV are presented in Figs. f|a] - [k| The nucleon decay 
amplitude is dominated by the wino dressing diagrams if r(wino) / r(total) w 1, while 
it is dominated by the gluino dressing diagrams if r (wino) /r (total) ^> 1. There 
occur cancellations between the wino dressing contributions and gluino dressing 
contributions in r(wino)/r(total) <C 1 region. The gluino contributions are small 

§This corresponds to taking Al ~ 0.46, where is the low energy QCD factor used in litera- 
tures H ||, H ^, It is argued the two-loop analysis gives A\, w 0.28 in Ref. jl6| 
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for any modes in the small tan f3 case. For large tan /3, however, gluino dominant 
region is realized in the charged lepton modes. This is brought about by the following 
two reasons: (1) in the simplified analyses based on an assumption of the diagonal 
gluino coupling, the gluino dressing diagrams contribute only to the Kv modes due 
to the color antisymmetry in the dimension-five coupling || |B| . In the present case, 
the gluino coupling is not diagonal any more and hence the gluino dressing processes 
contribute to the decay modes other than Kv modes. The nonvanishing off-diagonal 
gluino coupling in the up sector for large tan/3 significantly contributes to the 
charged lepton modes. (2) furthermore, if the squark masses are much larger than 
the gaugino masses, the amplitude from the gluino dressing diagrams is enhanced by 
a factor of (a 3 {mz) / ct^irnz)) 2 compared to that from the wino dressing diagrams, 
since F in Eq. ( [4.4| ) is asymptotically 

1 M 

F(m,m;M)~ - — for m > M , (5.4) 

167T 2 m 2 

and M 3 (mz) /M 2 (mz) = ct 3 {mz) I a^ijnz) because of the GUT relation of gaugino 
masses. Since ( |5.4|) gives an overall suppression of the nucleon decay amplitude for 
in ^> M (see Fig. |6^), the dominant gluino contribution in Figs. ^B] and |4c] is realized 
in the long lifetime region (see Fig. ^j. 

Translating the scanned parameters (mo, M g x, Ax) into the MSSM param- 
eters (rn^ L , M 2 , /i) where we take m^, the down-type squark mass of the first 
generation as the typical squark mass, we plot the calculated points in the MSSM 
parameter space for m top = 150 GeV and tan (3 = 2 in Figs. 5a and ^Fj. The region 



A in Fig. |6a| is excluded by LEP constraints on charginos and neutralinos pi 
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B(Z - 


>xx) 


< 


5 x 1(T 5 



where is a chargino, x is the lightest neutralino and x' is a heavier neutralino. 
No solution with radiative SU(2)xU(l) breaking is found in region B, which is a 
forbidden region. Points plotted with small dots are excluded due to the present 
lower bound for the proton lifetime r(p — * K + V) > 10 32 yrs, giving the constraint of 
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> 300 GeV. This constraint for [i is roughly unchanged for different m top and/or 
tan (3. Fig. |I5| shows the squark mass bound mg L > 400 GeV . If the lower bound 
for the proton lifetime is raised to r(p — > K + V) > 10 33 yrs with the near future 
experiment at Super-KAMIOKANDE, most of the parameter region with mg < 
1 TeV will be excluded. Lower bounds for other first and second generation squark 
masses are found similar to that for while the bound for the third generation 
squarks is lower in general due to the renormalization effect and the left-right mixing. 
Since the nucleon lifetime is approximately proportional to (tan/5) -2 0, the lower 
bound for the squark mass is raised for larger tan/3. In fact, r(p — » K + V) > 10 32 
yrs implies mj > 1 TeV for tan/3 = 30. 

VI Conclusion 

In this paper we have made a systematic analysis of the flavor mixing in the gaug- 
ino couplings within the framework of the minimal SUGRA-GUT. We have solved 
the one-loop RGEs for all MSSM parameters including off-diagonal Higgs coupling 
matrices with five input parameters, namely (m top , tan/3) at the electroweak scale 
m z and (m , M gX , A x ) at the GUT scale M x , and we have numerically obtained 
full mass spectra and mixing matrices, which satisfy the radiative electroweak sym- 
metry breaking condition. For a small tan/3 (tan/3 = 2), we have obtained a result 
consistent with the semi-analytic study [|TTJ, in which the top Yukawa coupling is 
assumed to be much larger than other Yukawa couplings: the left-left sector of 
the generation mixing matrix in the down-type quark-squark-gluino coupling is ap- 
proximately equal to the Kobayashi-Maskawa matrix, while the off-diagonal mixing 
matrix elements in the up-type gluino coupling are small. On the other hand, for 
large tan/3 = 10 and 30 where the bottom (and tau, for extremely large tan/3) 
Yukawa coupling is not negligibly small compared with the top Yukawa coupling, 
we have found that nonvanishing generation mixing in the up-type gluino coupling 
occurs with the magnitudes comparable to the corresponding Kobayashi-Maskawa 
matrix elements. The generation mixing in the down-type gluino coupling is also 
changed considerably. 

We have applied the generation mixing to the calculation of nucleon decay 
widths to study the contributions from the gluino dressing diagrams compared with 
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the wino dressing diagrams. In result, it is found that the gluino dressing diagrams 
give the dominant contribution to the decay mode containing a charged lepton if 
tan/3 ^> 1 and M3 <C rriq (typical squark mass). For the charged lepton modes 
with small tan /3, or the (anti-) neutrino emission modes with any tan (3, the gluino 
contributions are relatively small. In those cases, the contributions from the gluino 
dressing are at most of the same order of magnitude as the wino dressing contri- 
butions. We have scanned the MSSM parameter space to find allowed regions with 
the present constraints given by the nucleon decay experiments and the accelera- 
tor experiment s[]. The latter excludes the parameter region of small superpartner 
masses, and the former gives a strict bound to the masses of first and second gen- 
eration squarks. We argue that the whole parameter region with mj < 1 TeV in 
the minimal SU(5) SUGRA-GUT model can be tested by Super-KAMIOKANDE. 

Our method of calculations and the numerical result itself are adaptable to the 
analyses of FCNC in the minimal SUGRA model, which will be discussed elsewhere. 
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"in addition, cosmological constraint will be given by the analyses of the relic abundance of the 
lightest superprticle P2|. 
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Figure Captions 



FIG. |l| Examples of dressing diagrams which contribute to the proton decay process 
p — > K + V. 



FIG. [2a] Histograms of the mixing matrix element {Uu)\ (ul ~~ °l mixing) for 



tan = 2, 10 and 30 with m top = 150 GeV , 10 GeV < m < 10 TeV 
10 GeV < M gX < 10 TeV and -5 < A x < +5. 



FIG. 2B Histograms of the mixing matrix element {Uu)\ {ul — £l mixing) for 
tan/3 = 2, 10 and 30 with m top = 150 GeV . The parameters are same 



as those in Fig. 2a 



FIG. [2c| Histograms of the mixing matrix element (U' D )\ {di — sl mixing) for 
tan/3 = 2, 10 and 30 with m top = 150 GeV . The parameters are same 



as those in Fig. 2a 



FIG. 



Histograms of the mixing matrix element (U' D )\ [di — bi mixing) for 



tan/3 = 2, 10 and 30 with m t , 



op 



150 GeV . The parameters are same 



as those in Fig. |2a 



FIG. 



Histograms of the mixing matrix element {U' D )\ (sl — mixing) for 



tan/3 = 2, 10 and 30 with m t , 



op 



150 GeV . The parameters are same 



as those in Fig. |2a 



FIG. |2| Histograms of the mixing matrix element {U'd)\ {sl—^l mixing) for tan (5 



2, 10 and 30 with m t , 



op 



150 GeV . The parameters are same as those in 



Fig. 2a. 



FIG. |3| Nucleon decay partial lifetimes for m top = 150 GeV and tan/? = 2 with 
10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. 
The lepton in each mode has a lepton number —1 ('V, "e" and mean 17, 
e + and fi + , respectively). tl K" means a K meson with a s quark (K + or K°). 
V without a suffix means the total of three neutrinos. The shaded region is 
excluded experimentally. If the data points with r(n — > K°V) < 0.86 x 10 32 
yrs are omitted, the minimum value of each mode is raised to the vertical line. 
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FIG. |4a] The ratios of partial lifetimes calculated only the wino dressing diagrams 
(r(wino)) and those calculated with both wino and gluino dressing diagrams 
(r(total)) for m top = 150 GeV and tan/5 = 2 with 10 GeV < m < 10 TeV , 
10 GeV < M gX < 10 TeV and -5 < A x < +5. 



FIG. |4b] The ratios of partial lifetimes calculated only the wino dressing diagrams 
and those calculated with both wino and gluino dressing diagrams for m top = 
150 GeV and tan/3 = 10 with 10 GeV < m < 10 TeV , 10 GeV < M gX < 
10 TeV and -5 < A x < +5. 



FIG. |4c| The ratios of partial lifetimes calculated only the wino dressing diagrams 
and those calculated with both wino and gluino dressing diagrams for m top = 
150 GeV and tan/3 = 30 with 10 GeV < m < 10 TeV , 10 GeV < M gX < 
10 TeV and -5 < A x < +5. 

FIG. [5] A scatter plot of the p — > K + T> mode lifetime versus the wino-total ratio 
of the p — > K°e + mode for m top = 150 GeV and tan (3 = 10 with 10 GeV < 
m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. The shaded 
region is excluded experimentally. 



FIG. |6a] Scatter plots in /X-M2 plane for m top = 150 GeV and tan/3 = 2 with 
10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. 
Region A is excluded by LEP experiment. Region B has no radiative breaking 
solutions (see text). The region plotted with small dots should be excluded 
by the nucleon decay experiments. 



FIG. |6b] Scatter plots in mg L -M<i plane for m top = 150 GeV and tan/3 = 2 with 
10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. 
The region plotted with small dots should be excluded by the nucleon decay 
experiments. 
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Figures 



FIG. 1: Examples of dressing diagrams which contribute to the proton decay process 
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FIG. 2a: Histograms of the mixing matrix element (Uy)\ (ul — Cl mixing) for 
tan/? = 2, 10 and 30 with m top = 150 GeV , 10 GeV < m < 10 TeV , 10 GeV < 
M gX < 10 TeV and -5 < A x < +5. 



FIG. 2b: Histograms of the mixing matrix element {Uu)\ {ul ~ mixing) for 
tan (3 = 2, 10 and 30 with m top = 150 GeV . The parameters are same as those in 
Fig. EaL 



FIG. 2c: Histograms of the mixing matrix element (U' D )\ {di — mixing) for 
tan /3 — 2, 10 and 30 with m top = 150 GeV . The parameters are same as those in 
Fig. EaL 



FIG. 2d: Histograms of the mixing matrix element (6^)3 (di — bi mixing) for 
tan (3 = 2, 10 and 30 with m t0 p = 150 GeV . The parameters are same as those in 



Fig. £a 



FIG. 2e: Histograms of the mixing matrix element {JJ' D )\ (sl — di mixing) for 
tan (5 = 2, 10 and 30 with m top = 150 GeV . The parameters are same as those in 



Fig. ga 



FIG. 2f: Histograms of the mixing matrix element {0^)1 (s^ — b^ mixing) for tan/3 
2, 10 and 30 with m top = 150 GeV . The parameters are same as those in Fig. 
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FIG. 3: Nucleoli decay partial lifetimes for m t0 p = 150 GeV and tan/5 = 2 with 
10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. The 
lepton in each mode has a lepton number —1 ('V", "e" and "/i" mean 17, e + and 
respectively). li K" means a K meson with a s quark (K + or K°). V without a suffix 
means the total of three neutrinos. The shaded region is excluded experimentally. 
If the data points with r(n — > K°V) < 0.86 x 10 32 yrs are omitted, the minimum 
value of each mode is raised to the vertical line. 



FIG. 4a: The ratios of partial lifetimes calculated only the wino dressing dia- 
grams (r(wino)) and those calculated with both wino and gluino dressing diagrams 
(r(total)) for m top = 150 GeV and tan/3 = 2 with 10 GeV < m < 10 TeV , 
10 GeV < M gX < 10 TeV and -5 < A x < +5. 



FIG. 4b: The ratios of partial lifetimes calculated only the wino dressing diagrams 
and those calculated with both wino and gluino dressing diagrams for m top = 150 
GeV and tan/3 = 10 with 10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV 
and -5 < A x < +5. 



FIG. 4c: The ratios of partial lifetimes calculated only the wino dressing diagrams 
and those calculated with both wino and gluino dressing diagrams for m top = 150 
GeV and tan/3 = 30 with 10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV 
and -5 < A x < +5. 



FIG. 5: A scatter plot of the p — > K + v mode lifetime versus the wino-total ratio of 
the p — > K°e + mode for m top = 150 GeV and tan/? = 10 with 10 GeV < itlq < 
10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. The shaded region is 
excluded experimentally. 
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FIG. 6a: Scatter plots in 11-M2 plane for m top = 150 GeV and tan/3 = 2 with 
10 GeV < m < 10 TeV , 10 GeV < M gX < 10 TeV and -5 < A x < +5. Region 
A is excluded by LEP experiment. Region B has no radiative breaking solutions 
(see text). The region plotted with small dots should be excluded by the nucleon 
decay experiments. 



FIG. 6b: Scatter plots in -M 2 plane for m top = 150 GeV and tan/3 = 2 with 
10 GeV < m < 10 TeV , W GeV < M gX < 10 TeV and -5 < A x < +5. The 
region plotted with small dots should be excluded by the nucleon decay experiments. 
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